
function y = difWAD_fixedW(w,t,P,X,Z,R,C)

y = 0;
PZZ = calculatePHat_fixedW(w,t,X,Z,R,C);
for i = 1: t
  
    PP(:,:) = P(i,:,:);   
    PZZZ(:,:) = PZZ(i,:,:);
        
    for ii=1:R
        for jj=1:C
            y = y + myfunc(ii,jj,PP,PZZZ);
        end
    end
     
end

end



function v = myfunc(ii,jj,P,PZ)

p = P(ii,jj);
pHat = PZ(ii,jj);


v = abs(p-pHat)*pHat;

% if pHat ~=0
%     v = n * (p - pHat)^2 / pHat*(1-pHat) ;
% else
%     v = 0;
% end

end

